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Abstract 

A first principle reciprocating quantum refrigerator is investigated to determine the limitations of 
cooling to absolute zero. If the energy spectrum of the working medium possesses an uncontrollable 
gap, then there is a minimum achievable temperature above zero. Such a gap, combined with a 
negligible amount of noise, prevents adiabatic following during the expansion stage which is the 
necessary condition for reaching Tc — > 0. 
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Reciprocating refrigerators operate by a working medium shuttling heat from the cold to 
the hot reservoir. This requires external control of the temperature of the working medium. 
Upon contact with the cold side the working medium temperature has to become lower 
than Tc-the cold bath temperature. At very low temperatures a quantum description of the 
working medium is required where the control of temperature is governed by manipulating 
the energy levels of the system. 

A generic working medium possesses a Hamiltonian that is only partially controlled ex- 
ternally: 

H = H,„t + He,,H (1) 

where uj = uj(t) is the time dependent external control field. Typically, [Hj„(,Hext] 7^ 0., 
therefore, [H(t),H(t')] 7^ as a result a state diagonal in the temporary energy eigenstates 
cannot follow adiabatically. This fact, which is the source of quantum friction, has a profound 
effect on the performance of the heat pump l|, |2| . Almost perfect adiabaticity is the key to 
low temperature refrigeration. Typically, the internal interaction leads to an uncontrollable 
finite gap hJ in the energy level spectrum between the ground and first excited state. 
We will show that this gap combined with unavoidable quantum friction leads to a finite 
minimal temperature. An exception is a controllable energy gap which can be reduced to zero 
following the approach of the cold bath temperature Tc to zero. This case has been studied 
separately and leads to a vanishing rate of cooling when the absolute zero is approached . 

The Cycle of Operation, the Quantum Heat Pump. The working medium 
in the present study is composed of an interacting spin system. Eq. ([1]) is modeled 
by the SU(2) algebra of operators. We can realize the model by a system of two 
coupled spins tlmt = ^f^J {^l ® <t^ — (t^ ® o"^ ) = hJ'B2 where <t repre- 
sents the spin-Pauli operators, and J scales the strength of the inter particle interac- 
tion. For J — » 0, the system approaches a working medium with noninteracting atoms 
^. The external Hamiltonian represents interaction of spins with an external magnetic 
field: He^t = lhuj{t) (s-l 0p +P a^^ = cj(t)Bi. The SU(2) is closed with 
B3 = l{&l®&l + &l^al ) and [Bi,B2] = 2^63. 

The total Hamiltonian then becomes: 

H = (cu(t)Bi + JB2) . (2) 
The temporary energy levels, the eigenvalues of H are ei = —HQ, €2/3 = 0, €4 = HQ where 



Q = \^uj'^ + .P. For J 7^ there is a zero field splitting. Eq. ([2]) contains the essential 
features of the Hamiltonian of magnetic materials 5j. 

The dynamics of the quantum thermodynamical observables are described by completely 
positive maps within the formulation of quantum open systems . The dynamics is 

generated by the Liouville superoperator, £, studied in the Heisenberg picture, 

H = kP.au- £.(A) + ^ . (3) 

where £d is a generator of a completely positive Liouville super operator. 

The cycle studied is composed of two segments where the working medium is in contact 
with the cold/hot baths and the external control field uj is constant, termed isochores. In 
addition, there are two segments termed adiabats where the external field uj{t) varies and 
with it the energy level structure of the working medium. This cycle is a quantum analogue 



of the Otto cycle ^]. Each segment is characterized by a quantum propagator Us- The 
propagator maps the initial state of the working medium to the final state on the relevant 
segment. The four strokes of the cycle (see Fig. [1]) are: 

• Isochore (isomagnetic) A B\ the field is maintained constant uj = ujh the working 
medium is in contact with the hot bath of temperature T^. Cd leads to equilibrium 
with heat conductance F/j, for a period of t^. The segment dynamics is described by 
the propagator Uh- 

• Expansion adiabat (dimagnetization) B C: The field changes from Uh to uj^. in a 
time period of r^c- '^d = '^n represents external noise in the controls. The propagator 
becomes Uhc which is the main subject of study. 

• Isochore (isomagnetic) C ^ D: the field is maintained constant u = the working 
medium is in contact with the cold bath of temperature T^. Cn leads to equilibrium 
with heat conductance Fc, for a period of Tc- The segment dynamics is described by 
the propagator Uc- 

• Compression adiabat (magnetization) D ^ A: The field changes from Uc to Uh in a 
time period of Tch, jCd = jCn represents external noise in the controls. The propagator 
becomes Uch- 
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The product of the four propagators, Us is the cycle propagator: 



(4) 



Eventually, independent of initial condition, after a few cycles, the working medium will 
reach a limit cycle characterized as an invariant eigenvector of Ucyc with eigenvalue l(one) 
[2I. The characteristics of the refrigerator are therefore extracted from the limit cycle. 




FIG. 1: Refrigerator cycle in the frequency entropy plane. The von Neumann entropy Sm = 
—tr{plog p} (ABCD rectangle) as well as the energy entropy Se = —Y2Pi log Pi are shown {pi is 
the population of energy level i). The hot and cold isotherms are indicated. Both the expansion 
adiabat and the compression adiabats revolve exactly three periods. 

The dynamics of the expansion adiabat . The key to low temperatures is the expansion 
adiabat. In magnetic salt based refrigerators this segment is termed the adiabatic demagne- 



tization stage p , lid, lull • A necessary condition for cooling is that the energy of the working 
medium at contact point C is lower than the equilibrium energy at temperature Tc. What 
are the starting conditions at the beginning of the expansion segment point B? Consider- 
ing that the efficiency is limited by the Carnot cycle rjotto < Vcamot leads to the condition 
rp^. Now J and Qh > therefore using VLc{min) = J: 



(5) 
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This condition relates the hot end frequency Uh to the temperature Tc [5']. To force Tc to 
zero ujh ^ 00 and with it Qh- Under these conditions at equilibrium all population is in 



the ground state and {tl)h = —hQh- This is the optimal starting point for the expansion 
adiahat. On the cold side the necessary condition for refrigeration is that the internal energy 
of the working medium at the end of the expansion is smaller than the equilibrium energy 
with the cold bath. 

(H), < (H)e,(T,) = (l - 2e , (6) 

where (H)eg(Tc) is approximated by the low temperature limit ^ ksTc- Such a condition 
is fulfilled if the populations follows adiabatically the ground state during the expansion 
adiahat. Then (H)c = —Mlc and Eq. ([6]) is fulfilled. The expansion stage requires to reduce 
the external field uj from a large to a very small value maintaing adiabaticity. 

The orthogonal set of time independent operators Bi, is closed to the dynamics, and 
therefore they can supply a complete vector space to represent the propagators Uhc- A more 
thermodynamically oriented alternative is based on a time dependent set. The set includes 
the energy H and two other orthogonal operators: 



H = o;(t)Bi + JB2 , L = -JBi + uj{i)B2 , C = fi(t)B3 . (7) 

In general the dynamics on the expansion adiahat is generated hy C = Ch + where 
Ch = ■|[H, ■] and H()f:) the time dependent Hamiltonian Eq. ([2]). The external noise 
generator is Cn defined later. For perfect adiabatic following the propagator U factorizes 
between H, and L and C. 

The noiseless dynamics generated only by the Hamiltonian H(t) is the key to adiabaticity: 
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The ability of the working medium to follow the energy spectrum is defined by the adiabatic 
measure /i = If /x = the propagator factorizes. Constant /i minimizes the non- 

adiabatic deviations during the expansion. In addition constant yU leads to a closed form 
solution for the propagator Uhc forcing a particular scheduling of the external field loit) with 



+ 



time: ujit) = J f / a/1 — where/ is a linear function of time: f — ^ 
The adiabatic parameter /i and the time allocated to the adiahat The, obey the reciprocal 



relation: /x = where Khc 
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Eq. ([H]) is integrated by defining a new time variable: dO = Qdt The final values of Qhc 
becomes: Qhc = Thcj^^hc where: ^hc = ^arcsin(^) — arcsin(^) j and > $ > — |. 

Eq. (IHl) is solved by noticing that the diagonal is a unit matrix multiplied by a time 
dependent scalar. Therefore we seek a solution of the type Uhc = U1U2 where \Ui,lA2\ = 0. 
The integral of the diagonal part of Eq. ([8] ) becomes: 

= e^Io'^^^'^I = , (9) 

which can be interpreted as the scaling of the energy levels with the variation in Q. 

To integrate U2 the non diagonal parts of Eq. ([H]), are diagonalized, leading to the 
eigenvalues 0, —iy^, i^-, where q = a/1 + /x^, and the propagator: 



El c 



At(l-c) s M^+c 
\ g2 q j 



(10) 



where s = sin{qQ) and c = cos{q<d). 

The adiabatic limit is described by /x — 0. Then Eq. ( ITOi) factorizes. These are the 
perfect adiabatic following conditions. In general Eq. ffTOj) describes a periodic motion of H 
L and C. Each period is defined by 

qQ = 27il 1 = 0,1,2... (11) 

where / is the winding number. At the end of each period U2 restores to the identity matrix. 
These are the frictionless conditions of adiabatic following. For intermediate times (H) is 
always larger than the frictionless value. The amplitude of this periodic dynamics decreases 
when m becomes smaller, Cf. ^2(1, 1) in Eq. ( ITOi) . 

The frictionless conditions define a quantization condition for the adiabatic parameter /x: 

Examining Eq. (fT2l) we find that there is no solution for / = 0. The first frictionless solution 
I > V ^ leads to a minimum expansion time for frictionless solutions: 




Thc{min) = Khc\ { ^] - 1 • (13) 



The family of all frictionless solutions leads to refrigeration cycles which obey Eq. for 
any > 0. Such frictionless refrigerators have no minimum temperature above = 0. 

The Effective Minimal Temperature . Any realistic refrigerator is subject to noise 
on the external controls. Perfect adiabaticity requires precise control of the scheduling of the 
external field uj{t). Any deviation from perfect adiabatic following, maintaining the ground 



state on the expansion adiabat will lead to a minimum temperature. If (H) 
where 6 is the deviation from perfect adiabatic following then, from Eq. 
leading to: 

Tr > . / . > 



6 < 2e 



hJ 



(14) 



h\og{5l2) - -A;,log((5/2) ' 
We will now show that even an insignificant amount of noise will lead to Tc{min) > 0. 

First we consider a piecewise process controlling the scheduling of u in time. At every time 

interval, u is updated to its new value. Then random errors are expected in the duration of 

these time intervals described by the Liouville operator C^. This process is mathematically 

equivalent to a dephasing process on the expansion adiabat [12 1. This stochastic dynamics 

can be modeled by a Gaussian semigroup with the generator [8, 



-|[H,[H,A] 



(15) 
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which is termed phase noise. The modified equations of motion on the adiabats become: 

/w\ /a n 

- 7pfi -1 

1 §-ipn^ 

We seek a product form solution: Uhc = U1IA2IA3 where W3 the noise propagator, given by: 
^U,{t) = W{t)Us{t), where: 
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W{t) = U2{-t) 
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(17) 
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We seek an approximate solution for W3 in the limit when /i — > 0, then U2 = 1 since this is 
the frictionless limit. Expanding Eq. ( |T71) to first order in /i leads to: 
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UsiThc) is solved in two steps. First evaluating the propagator for one period of 6: for which 
Q{t) is almost constant, and then the global propagator becomes the product of the one 
period propagators for / periods: U-s^Thc) ~ ^3(6 = 27r)'. The Magnus expansion to second 
order is employed to obtain the one period propagator ^3(27?): 

W3(e = 27r) ^e^^+^^+- (19) 

where: Mi = J^"" dQW{e) and M2 = ^ /q^" jf dQdQ'[W{e), W{e')] + ...)• The first order 
Magnus term leads to: 

/ 1 

^3(6 = 27r)M, ^ 

which to first order in /i, 6 the deviation from perfect adiabatic following is zero. The second 
order Magnus approximation leads to: 



/i(i - e-2^>^; 

g-27r7pC 



(20) 
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where S = sin a and C = cos a. a = iT'^pVl^^J^^'^ + 4 and as ;U ^ 0, a = 27r7pi7/i. The 
second order propagator Us{Thc), for I revolutions is also a rotation matrix identical to Eq. 



(I2TI) . with a new angle ai = 27r7pyU J^'''' ^l{Q)dQ = 717^ J In 
of U3 from the identity operator defines 6. Asymptotically for /i ^ , 6min 
7r^7p ln[co'/i/ J]. Any time variation in fi will lead to 5 > 6min- 

Another source of external noise is due to amplitude errors in the control of the frequency 
uj{t). These errors are modeled by a Gaussian random process described by the Lindblad 
term: C]\f{A) = — 7aCo'^[Bi, [Bi, A]] where 7a characterizes the amplitude noise. The equa- 
tion of motion for the noise propagator U3 becomes: ^W3(6) = W(0) ^3(6) where: 
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We seek an approximate for small /i. Using a similar procedure of calculating the propagator 



for one period ^3(0 = 27r). The ^3(1, 1) element decouples from the remaining part of the 



propagator. As a result 1 — 6 = ^3(1, 1) ~ e "'^''^''''"1. The smallest 6 is achieved for a one 
period cycle, Eq. ([T3l) then 5min ~ 47aJT^. 

Figure [2] shows 5 as a function of propagation time for different values of /i calculated 
numerically for phase (integrating Eq. ffTSl) ) and amplitude noise (integrating Eq. ([22]) ). 




FIG. 2: The deviation 5 from perfect adiabatic following as a function of the expansion time for 
different revolutions / for a) phase noise and b) amplitude noise. The end points of each expansion 
are marked by a filled circle corrosponding to the value of 6 at the end of the expansion adiabat. 
For the phase noise the final S decreses with The while for amplitude noise 5 increases with The- 

Figure [S] shows the minimum temperature calculated numerically as a function of ex- 
pansion time Tfic, for the phase and amplitude noise. The exact numerical calculation are 
consistent with the approximation when fi —>■ 0. The phase noise has a monotonic decrease 

of TJmin) reaching saturation as The 00 where TJmin) = 7^ -^7. TJmin) 

of the amplitude noise is monotonically increasing function of time which means that short 
expansion times lead to the minimum temperature. If both amplitude and phase noise op- 
erate simultaneously the minimum temperature will be obtained at the crossing point. This 
optimum will move with the ratio 7p/7a- 
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FIG. 3: The minimum temperature as a function of the expansion adiabat time ahocation The 
(bottom scale) and / the winding number (upper scale), for the two noise models. 0,fi = 3.8, 

-0.001 and 7a = -0.0001. 



= 2.01, J = 2 7p 



Conclusions: The necessary condition for the working medium to cool down to absolute 
zero is that (H)c = —hQc- In order to start in the ground state at the hot end the working 
medium has to equilibrate with a very high frequency Qh- Perfect adiabatic following will 
maintain the system in its ground state which defines the frictionless solution. We found 
a family of additional frictionless solution obeying a quantization rule for /i the adiabatic 
parameter. 

The main result of this study is that any noise in the controls of uj{t), will eliminate the 
frictionless solutions leading to a minimum temperature Tc{min). This finding is consistent 
with experiments on demagnetization cooling of a gas 13] which obtained a minimum tem- 



perature an order of magnitude larger than the theoretical prediction 15| attributing the 
discrepancy to noise in the controls. The logarithmic dependence on the noise parameters 
means that Tc(mm) is of the oder of J. The surprise is the negative effect of phase noise. 
Its generator Eq. ( fTSi) can be interpreted as the result of weak measurement of H, forcing 
collapse to a state diagonal in energy. It would seem that this should improve adiabaticity 
and reduce friction. For an engine model we found such an effect of enhanced performance 



12l |. Nevertheless even when /i is small and constant, which is the condition for minimal 
phase noise, the non-commutativity of H leads to a finite 6min which saturates as r^c ^ oo, 
and therefore TJmin) > 0. 
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